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Owing to its remarkable electronic and transport properties, graphene has great potential of replacing 
silicon for next-generation electronics and optoelectronics; but its zero bandgap associated with Dirac 
fermions prevents such applications. Among numerous attempts to create semiconducting graphene, 
periodic patterning using defects, passivation, doping, nanoscale perforation, etc., is particularly promising 
and has been realized experimentally. However, despite extensive theoretical investigations, the precise role 
of periodic modulations on electronic structures of graphene remains elusive. Here we employ both the 
tight-binding modeling and first-principles electronic structure calculations to show that the appearance of 
bandgap in patterned graphene has a geometric symmetry origin. Thus the analytic rule of gap-opening by 
patterning graphene is derived, which indicates that if a modified graphene is a semiconductor, its two 
corresponding carbon nanotubes, whose chiral vectors equal graphene's supercell lattice vectors, are both 
semimetals. 

Graphene has risen as a fascinating system in condensed matter physics for not only fundamental science 
but also technological applications 1 4 because of its many unique and amazing properties originating 
from its massless fermions 2,5 due to the linear dispersion of energy bands near the Dirac points (K). For 
example, under gate voltage charge carriers in graphene can be tuned continuously between electrons and holes 1 , 
and their mobilities exceed 2.5 X 10 s cm 2 V _1 s _1 under ambient conditions 6 . This certainly makes graphene 
extremely attractive for electronic and optoelectronic applications 3,4 ; however, graphene is a semimetal with 
bandgap (£ g ) closing at K, which prohibits switching off the graphene channels in field-effect transistors and 
building functional junctions in graphene optoelectronics. 

Extraordinary efforts 716 have been spent to create semiconducting graphene materials that maintain the 
exceptional transport property. One scheme is obtaining bandgap through quantum confinement in graphene 
nanoribbons 7 (GNRs). But it is difficult to control the magnitude of bandgap in GNRs, since £ g is very sensitive to 
ribbon edge and width; moreover, a narrow GNR, which is necessary for a sizable gap, cannot carry sufficiently 
large currents. Other methods include substrate-induced 8,9 and strain-induced 10 bandgaps, but the former is 
essentially not tunable, while the latter is practically limited by the ability of tuning £ g to the desired value. An 
alternative approach is to periodically modify graphene 17 , which is crucial to preserve its anomalously high charge 
carriers mobilities, by normal or inverse Stone-Thrower- Wales type of defects 11,12 , hydrogen passivation 13 , boron 
and nitrogen doping 14 , and nanoscale holes creating graphene nanomeshes 1517 (GNMs), etc. In particular, Balog 
et al." demonstrated the existence of bandgap in graphene by patterned hydrogen adsorption, and Bai et al. 16 
made functional GNM field-effect transistors capable of supporting currents 100 times of those in GNR-based 
devices. 

However, it is still a mystery that a periodically modulated graphene can either open up a substantial gap or 
remain gapless (or a tiny gap) if its supercell lattice changes only slightly 18 28 , which complicates the creation of 
bandgap by patterning and is apparently at odds with the proposal that bandgap in these graphene structures is 
due to quantum confinement between neighboring modulated sites 16 . The quantitative description of quantum 
confinement that results in energy gap in a patterned graphene is missing; instead, previous theoretical investi- 
gations 1828 mainly focused on some special patterns whose supercell lattices form rectangular or hexagonal 
structures. A number of interesting empirical rules have been unraveled, but a fundamental and precise under- 
standing on how electronic structures are affected by regularly patterning graphene has not yet been revealed. 

Results 

Analytical theory. In this work we adopt and extend the tight-binding model 2 , which describes the sp 2 -bonding in 
graphene extremely well, to include periodic structural alterations modeled by an external potential I7(r) applied 
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on graphene 29 . We denote the primitive lattice vectors of the pristine 
graphene by and a 2 , and the supercell lattice vectors of the 
patterned graphene can be written as 

Ri = «]ai + miH 2 

R 2 = n 2 ai + m 2 a 2 

with (n 1 , nti, n 2 , m 2 ) four integers. If the perturbed sites form two- 
dimensional Bravais lattices, then the external potential (7(r) has the 
same translational symmetry 

U(t) = U(r+l l R 1 + kR2), (2) 

where l x and l 2 are two integers. 

At Dirac points K, the intrinsic tight-binding off-diagonal matrix 
elements /.(/(k) of the pristine graphene vanishes, giving rise to a zero 
bandgap; therefore, a non-vanishing U(K), the Fourier transforma- 
tion of U(r) at Dirac points will result in an energy gap. This leads to 
the following bandgap opening rule: 

J «i — mi=3p and n 2 — m 2 = 3q, (Gap opening) 
l_ Otherwise, (Gap closure) 

where p and q are two integers. Our modeling is briefly summarized 
in the Methods section, and the physical reasoning and mathematical 
derivations are detailed in Supplementary Information. 

Equation (3) indicates that only one ninth of arbitrarily patterned 
graphenes are semiconducting while the rest remain semi-metallic; 
thus the patterning periodicity must be carefully chosen to create a 
gap. Furthermore, Equation (1) suggests that a patterned graphene 



can be mapped to a pair of carbon nanotubes 30 (CNTs) with chiral 
vectors Q, equal to Ri and R 2 , respectively; then bandgap opening in 
a patterned graphene has an exactly opposite correlation to that in 
CNTs: it possesses a sizable bandgap only if both of its corresponding 
CNTs are gapless at the tight-binding level. 

We illustrate our analytical theory by discussing a crucial and 
representative class of patterned graphene, namely GNMs. Fig. 1 
plots various unit cells of GNMs, which are parallelograms with 
different angles, side lengths, and orientations with respect to the 
pristine graphene lattice. Because of the honeycomb structure of 
graphene, 60° and 90° (rectangle) parallelograms are especially 
important and have also been realized experimentally 15,16 . Like other 
patterned graphenes, a GNM can be mapped to two CNTs. Since 
unfolding a CNT generates a GNR, a GNM can also be mapped to a 
pair of GNRs, as illustrated in Figs, lg and lh. Experimentally, GNRs 
have been created by unzipping CNTs longitudinally 31 , whereas 
CNTs could also be produced by wrapping GNRs around 32 . A 
GNR and its corresponding CNT have similar dependency of elec- 
tronic band structure on their transverse vectors, and the well-known 
three-fold alternation in electronic structures is not surprising, given 
the hexagonal symmetry of graphene 2 . 

However, a GNM (or other types of patterned graphene) 's 
dependence of bandgap opening on supercell lattice vectors is oppo- 
sitely to that of its corresponding CNTs (or GNRs), because along 
each lattice vector of a GNM the matching condition on Dirac points 
is exactly the same as the circular confinement condition of a CNT. A 
CNT is characterized by its chiral vector defined by a pair of integers 
(«, m): Q, = wai + ma 2 . The tight-binding model asserts that a CNT 
is semi-metallic if n — m = 0 mod 3, or semiconducting otherwise 30 , 




Figure 1 | Unit cells of GNMs. (a), (b) and (c) 60° parallelograms, with two lattice vectors along the zigzag (blue), armchair (gray) and chiral (green) 
directions, respectively, (d) Oblique parallelogram (red), (e) and (f) Rectangles whose lattice vectors are along the zigzag and armchair directions, 
respectively (yellow), or along two chiral directions (purple). Dangling bonds in GNMs and GNRs are passivated by hydrogen atoms (small gray dots), 
while carbon atoms are denoted by brown dots, (g) and (h) The relationship among GNMs, GNRs and CNTs, as explained in the text. 
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since its allowed k points satisfy k • Q, = 2nl with I integers. 
Figs. 2d-f illustrate that only when these lines of k points pass 
through graphene's Dirac points K, a CNT remains semimetallic. 
On the other hand, as shown in Figs. 2g-i, the reciprocal lattice 
vectors of a GNM are the intersections of two sets of parallel lines 
representing the allowed k points of its two matching CNTs. To 
create a bandgap, the reciprocal lattice vectors of a GNM must 
overlap the Dirac points, requiring both sets of parallel lines pass 
through K; consequently, a semiconducting GNM corresponds to 
two semimetallic CNTs. 

First-principles calculations. In the following we demonstrate that 
the above conclusion drawn from analytical modeling is in excellent 
agreement with first-principles results. We present electronic struc- 
tures obtained using the density functional theory (DFT 33 ' 34 ) for 
GNMs; in addition, our DFT calculations on periodically H-passi- 
vated and boron-nitrogen doped graphene show similar results, 
which will be published elsewhere. Although it is well known that 
DFT can severely underestimate bandgaps in semiconductors and 
insulators, DFT makes good predictions on the trends of electronic 



band structures and bandgap 35 . In particular, DFT correctly predicts 
a zero bandgap (closing at the Dirac points) for pristine graphene 2 , 
and the predicted trends for bandgaps and band structures in 
graphene-based nanostructures, such as CNTs and GNRs, are 
mostly correct, compared with the many-body perturbation theory 
within the GW approximation and experimental data 7,30,36 . If these 
materials are semimetallic (or nearly semimetallic), DFT predicts 
them to have a zero (or nearly zero) bandgap; on the other hand, if 
they are semiconductors, DFT predicts sizable £ g , which tend to be 
underestimated. 

Figs. 2j-l compare DFT electronic band structures of three oblique 
GNMs with (ni, mi, n 2 , m 2 ) = (9, -3, -4, 9), (7, 4, -4, 8) and (9, -2, 
— 4, 9), whose crystal structures are shown in Figs. 2g-i, respectively. 
Although these three GNMs have same holes and similar unit cells, 
only the middle one has a substantial gap because its two matching 
CNTs are semimetallic. The left GNM corresponding to one metallic 
and one semiconducting CNT, has a tiny gap at M 2 and it is nearly 
semimetallic; the GNM on the right remains gapless, with both of its 
matching CNTs semiconducting. These first-principles results agree 
well with our analytic theory, clearly demonstrating the sensitivity of 




Figure 2 | Band gap opening/closure in GNMs and their corresponding CNTs. (a), (b) and (c) Crystal structures of a pair of CNTs in each panel, 
whose chiral vectors are equal to the lattice vectors of GNMs plotted on the left side of (g), (h) and (i), and whose allowed k-points (parallel lines) are 
plotted in (d), (e) and (f), respectively. The right side of (g), (h) and (i) show the reciprocal lattices of GNMs plotted in (a), (b) and (c), respectively. In 
(d)-(i) the hexagons indicate the Brillouin zone of the pristine graphene, with 6 corners the Dirac points, (j), (k) and (1) Electronic band structures of 
GNMs plotted in (g), (h) and (i), respectively. Brillouin zones of all GNMs in the text are plotted in Supplementary Fig. 2. 
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band structure to patterning periodicity, whose small variation could 
modify electronic structures significantly. They also suggest that in a 
GNM the top valence bands and bottom conduction bands derived 
from the n orbitals are nearly symmetric with respect to the Fermi 
energy (or the middle of energy gap), consistent with the tight-bind- 
ing model for graphene in which only the nearest-neighbor interac- 
tions are taken into account. 

Next we discuss rectangular and hexagonal GNMs, which have a 
one-to-one correlation to CNTs, because their lattice vectors R l and 
R 2 are not independent. R : and R 2 of a rectangular GNM form the 
chiral and translational (T) vectors of a CNT subject to an integer 
multiplier j: Ri = and R 2 = ft, or Ri = ft and R 2 = Q,. Here T = 
t t ai + f 2 a 2 , with ( 1 = (2m + «)/(,t 2 = — (2n + m)/f, and t the greatest 
common divisor of 2m + n and 2n + m. If n — m is not a multiple of 
3, then t l — t 2 = 0 mod 3; thus at least n l — m 1 or n 2 — m 2 is divisible 
by 3. We set n 2 — m 2 = 0 mod 3, so that the existence of a bandgap in 
a rectangular GNM depends only on R t by the difference between Wj 
and nil. it is semiconducting if n t — m t = 0 mod 3 or semi-metallic 
otherwise, exactly in contrast to a CNT with C h = R it as summarized 
in Fig. 3i. The DFT band structures of two rectangular GNMs 
(Figs. 3a and 3b) with (ni, mi, n 2 , m 2 ) = (6, 3, —4, 5) and (6, 4, 
— 7, 8) are plotted in Figs. 3e and 3f. In both cases n 2 — m 2 are 
divisible by 3; the former is a semiconductor with a direct gap at 
the r-point, while the latter is semimetallic with gap closing at the 
T—M 2 segment. In fact, all semiconducting GNMs open up a direct 
bandgaps at the F point. 

In a rectangular GNM, ^ and R 2 are perpendicular to each other 
with conjugate chiralities so that at least one of its two matching 
CNTs are metallic, and then one third of rectangular GNMs are 




semiconducting. On the other hand, a 60° parallelogram GNM's 
lattice vectors R : and R 2 have the same chirality subject to integer 
multipliers. In general we cannot assign a pair of integers to char- 
acterize a 60° parallelogram GNM, because Ri and R 2 might be 
different in length. For Ri and R 2 with the same length, a 60° par- 
allelogram becomes a 60° rhombus and such GNMs have the hexa- 
gonal symmetry, in which n l — m 1 and n 2 — m 2 are either both 
divisible or neither divisible by 3. Therefore the existence of a band- 
gap in a hexagonal GNM also depends only on R l by n x — m lt and 
one third of them possess bandgaps, compared with only one ninth 
of 60° parallelogram GNMs semiconducting. Figs. 3g and 3h plot 
DFT band structures of two hexagonal GNMs (Figs. 3c and 3d) with 
(«!, mi, n 2 , m 2 ) = (5, 4, 9, —5) and (7, 1, 8, —7), respectively. As 
predicted by Eq. (3), the former is semimetallic, while the latter is a 
semiconductor. 

Previous theory and simulations 18-28 studied many special GNMs, 
namely the zigzag (Fig. la), armchair (Fig. lb), and zigzag- armchair 
(Fig. le) GNMs, in which R 1 and R 2 are both along the zigzag or 
armchair direction, or along zigzag and armchair directions respect- 
ively. A zigzag GNM is characterized by a pair of integers (Z x , Z 2 ) 
with mi, n 2 , m 2 ) = (Z lt 0, 0,Z 2 ), an armchair by (A lt A 2 ) with (n 1( 
mi, n 2 , m 2 ) = (2Ai, — Ai, A 2 , A 2 ) and a zigzag-armchair by (Z, A) 
with (rii, mi, n 2 , m 2 ) = (Z, —Z, A, A). The bandgap opening rules for 
these special cases as functions of structural parameters A and Z are 
presented in Figs. 4e, 4f and 4j, respectively. One third of zigzag- 
armchair GNMs are semiconducting, which is independent of the 
value of A and depends only on whether Z = 0 mod 3. A zigzag GNM 
acquires a bandgap if both Z l = 0 mod 3 and Z 2 = 0 mod 3, so only 
one ninth of them are semiconducting, whereas all armchair GNMs 




/7j — 6, m, — 3, rL, — — 4, m 1 — 5 



f «, = 6,m,=4,«,=-7,/n,=8 9 «, = 5, m| =4.n 2 =9.m : =-5 h ", = 7, m, = I.;,, = 8, = -7 




Figure 3 | Rectangular and hexagonal GNMs. (a) and (b) Crystal structures of two rectangular GNMs; (e) and (f) the corresponding electronic band 
structures, (c) and (d) Crystal structures of two hexagonal GNMs; (g) and (h) the corresponding electronic band structures, (i) Bandgap opening (red) 
and closure (blue) for rectangular and hexagonal GNMs (hexagons) and CNTs (circles) characterized by (n lt mi) for vector R l of GNMs and C;, of CNTs, 
respectively, i.e., C;, = R| = n^i + m^. Here the solid arrows in (i) denote the zigzag direction, while the dashed line indicates the armchair vector. 
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are semiconductors regardless of the values of A x and A 2 . These 
predictions have been verified by previous 18-28 and present DFT cal- 
culations whose results are plotted in Figs. 4a-d, 4g and 4h. 
Comparing band structures plotted in Figs. 2-4, one finds that for 
the same symmetry of patterning, semiconducting graphene alter- 
nates with semimetallic ones, except for the armchair case where all 
of them are semiconductors. Apparently, the armchair structure is 
the preferred choice to open bandgap by patterning graphene. 

We have carried out additional DFT calculations for electronic 
structures of other patterned graphenes with a variety of unit cells 
and defect sizes other than what reported in Figs. 2-4. All these 
results (in Supplementary Information) agree with Equation (3) on 
bandgap existence. Furthermore, the magnitude of bandgap of a 
semiconducting patterned graphene depends on the onsite perturba- 
tion strength and inter-site distances; in particular, £ g of GNMs as 
functions of hole size and unit-cell area can be well fitted to a simple 
scaling rule proposed by Pedersen et al. ls , as demonstrated in Fig. 4i. 
Hence once a semiconducting patterning is chosen, £ g can be tuned 
to the ideal values for target applications by adjusting the types and 
sizes of and distances between structural modulations. 

Discussion 

Using a perturbative tight-binding model we have derived a universal 
bandgap opening rule for graphene with arbitrary Bravais periodic 



patterning. A simple reverse relation on existence of bandgap 
between a modified graphene and its two analogous CNTs is found, 
which is confirmed by first-principles electronic structure calcula- 
tions. Current modeling can be easily generalized from Bravais to 
non- Bravais supercells, such as the honeycomb structure, and thus 
the DFT results for the honeycomb GNMs 26 can also be well 
explained. We conclude that bandgap opening in patterned graphene 
is not directly caused by quantum confinement as in CNTs or GNRs, 
though confinement still plays a role in determining the magnitude 
of £ g once a gap is present; instead, it originates from geometric 
symmetry breaking at the Dirac points, leading to the exactly oppos- 
ite bandgap opening conditions between modified graphene and 
CNTs. 

Thus our present work has solved an outstanding problem on how 
to properly pattern graphene to induce bandgaps, pointing to prac- 
tical routes of geometrically designing graphene with controllable 
and desirable electronic behaviors. Previous experiments 1317 have 
demonstrated that patterning graphene could open up a substantial 
bandgap suitable for transistors, but none of them have yet reported 
the predicted sensitivity of electronic behavior to patterning peri- 
odicity. Their patterned graphene structures are rather rough, which 
might greatly hamper transport of charge carriers, while a high-qual- 
ity accurate periodic patterning is expected to largely maintain the 
ultrahigh electron mobility. This requires the ability of tailoring 




A=7,Z=12 
Z,=12,Z 2 =12 
A.=7,A. 



Figure 4 | Zigzag, armchair, and zigzag-armchair GNMs. Electronic band structures of two zigzag (a, b), armchair (c, d), and zigzag- armchair 
(g, h), GNMs. Their crystal structures are plotted in Supplementary Fig. 3. (e), (f) and (j). Band gap opening (red) and closure (blue) for zigzag (e), 
armchair (f) and zigzag-armchair (j) GNMs characterized by [Z x , Z 2 ), (A x , A 2 ) and (A, Z), respectively, (i) Bandgap (eV) as a function of hole size 
indicated by the number of removed carbon atoms (AT relnove< j) per unit cell for a zigzag (red), an armchair (green) or a zigzag-armchair (black) GNM. 
Computed data (closed circles) are fitted to smooth curves, as described in Supplementary Information. 
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graphene structures precisely at the atomic scale 37 - an exciting and 
rapidly advancing frontier in nanoscience and nanotechnologies for 
realizing electronics and optoelectronics towards ever smaller length 
scales. 

Finally, the restriction of our modeling, as well as the key to its 
success, is that the valence and conduction bands in graphene touch 
only at the Dirac points. In other materials with vanishingly small 
density of states near the Fermi energy or, equivalently, bandgap 
closing only at a small discrete set of k points, it is possible to perturb 
all these touching k points by patterning. Moreover, it is appealing to 
investigate the effects of periodic perturbation and structural sym- 
metry breaking on strongly-correlated materials showing intriguing 
quantum/topological orders, such as topological insulators 38 , whose 
gapless surface states are protected by time-reversal symmetry and 
the topology of the material. 

Methods 

First-principles calculations. Our first-principles electronic structure calculations 
based on DFT were performed employing the SIESTA package 39 using a triple-^ 
polarized atomic basis set for carbon atoms, whose numerical accuracy has been 
rigorously examined against the planewave-based VASP program 40 . The two- 
dimensional Brillouin zones of GNMs were sampled on a 2 X 2 Monkhorst-Pack k- 
grid with good convergence due to very large unit cells. The generalized gradient 
approximation 41 was used for the exchange-correlation functional. All calculations 
are spin-polarized. GNM structures are optimized until the maximum atomic force is 
less than 0.02 eV/A. Dangling bonds in GNMs are passivated by hydrogen atoms. 

Tight -binding model. Graphene has two sublattices (A and B) whose Hamiltonian is 
written as 



H K = 



where E Q and x 0 are two parameters. The eigenvalues are 

£ ± (k) =£„+;.„ [/(k)|, 



^3 + 2cos[-/3k y aJ - 



■4 cos 



— * 
2 



-k,a 



(4) 



(5) 



(6) 



with a «* 1.42 A, the neighboring carbon-carbon distance. The structural 
modifications are modeled by an external periodic potential U(r). If the scattering 
matrix element between two degenerate pseudo-spin states at the Dirac points 



;. 1 (K)=N{^|t/(K)|^>, 



(7) 



is non-zero, where <j> A and <j) B are orbitals for a pair of the nearest neighbors, N is the 
effective average number of the nearest neighbors per site, and [7(k) is the Fourier 
transformation of U(r), or the scattering matrix element between two degenerate 
states at Dirac points K and K' 42 - 43 



<K| V, |K'> | 2 = | U(K - K) | 2 cos 2 (five/2) , 



(8) 



is non-zero, where 9 KK ' is tne angle between K and K', then there is a sizable bandgap 
in the defected graphene. Eqs. 7 and 8 lead to the same bandgap opening condition of 
C7(K) ^ 0. 
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